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Abstract 

On 

We present a lattice computation of the effective potential for 0(2)-invariant 
Q_i' (A<3? 4 )4 theory in the region of bare parameters corresponding to a classi- 

cally scale-invariant theory. As expected from "triviality" and as in the one- 
^ ; component theory, we find very good agreement with the onedoop prediction, 

while a perturbative leadingdog improvement of the effective potential fails 
to reproduce the Monte Carlo data. The mass rrih of the free shifted ra- 
dial field is related to the renormalized vacuum expectation value vr through 
the same relation = 87r 2 u|j as in the one-component case. This confirms 
the prediction of a weakly interacting 2.2 TeV Higgs particle in the standard 
model. 

In the study of Spontaneous Symmetry Breaking (SSB) the simplest quantity to compute 
is the vacuum expectation value of the bare scalar field $b(x) in the presence of an external 
constant source J 
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(*b)j = MJ) (1) 
Computing ^ on the lattice at several J- values is equivalent [IJ to inverting the relation 

J = J^b) = ^ (2) 



involving the effective potential V^//(0b). In this framework, the occurrence of SSB is 
determined from exploring (for J ^ 0) the properties of the function 

MJ) = -M-J) (3) 

in connection with the limiting behaviour at zero external source 

jim MJ) = ±%/0 (4) 
over a suitable range of the bare parameters (r G , A D ) in the lattice action 



A lattice simulation of the weakly coupled one-component massless A$ 4 theory has shown 
that there is a well defined region in the bare parameter plane (r Q , A D ), corresponding to SSB 
from the classically scale-invariant case, where the effective potential is reproduced by its 
one- loop form to very high accuracy. Even though ^§ << 1, a "leading- log improvement" 
of the one-loop potential completely fails to reproduce the Monte Carlo data. This result, 
while contradicting the naive perturbative expectations, confirms the crucial insight of ]3],|4| 
on the basis of the generally accepted "triviality" of (A$ 4 )4: for a "trivial" theory the 
one-loop potential is effectively exact in the continuum limit. In fact, since there are no 
observable interactions, the effective potential is just the sum of the classical potential and 
the zero-point energy of the free fluctuation field h(x) = &b{ x ) — (<&b)- The traditional per- 
turbative renormalization, based on the concept of a cutoff-independent and non-vanishing 
renormalized coupling at non-zero external momenta Xr, is not appropriate just because 
it would spoil this exactness-it does not properly re-absorb infinities but merely pushes them 
into "higher orders" which are then neglected. In this sense, the unphysical features of the 



perturbative renormalization (the one-loop Landau pole or the spurious two-loop ultraviolet 
fixed point at non zero bare coupling which contradicts the rigorous arguments of |J) are 
just a signal of the inadequacy of the procedure. However, it is simple to renormalize the 
one-loop potential exactly PH|7]-|TTf (this was first discovered in the context of the Gaussian 
effective potential JT^j-Tnj]). The particle mass m h is related to the ultraviolet cutoff and to 
the bare vacuum field vb through 

ry ry ry . 1 6 7T . / \ 

m h = 3X v B = A exp(— — ). (6) 

[Note that A Q = A^/6 in the notation of H§|]. At the same time, the vacuum energy density 
(a renormalization-group-invariant quantity) is 

V^°°*(±v B ) = (7) 

Hence, to get a cutoff-independent rrih in the continuum limit A — > oo, X Q = A G (A) has to 
vanish as l/(ln— ) and the bare vacuum field 4>b is non-trivially rescaled with respect to 
the physical field 4>r through 

01 = Z^ 2 R (8) 



with Zcj, ~ 1/A Q . The non-perturbative nature of the vacuum field renormalization 
(Z^ ~ 1/A G ), first discovered in the gaussian approximation by Stevenson and Tarrach 
||13|| , should not be confused with the /i-field wave function renormalization. Since h is just 
a free field, one has trivially Z^—l. The structure Z^ ^ Zh is allowed because for a scalar 
field the decomposition into = and ^ components is Lorentz-invariant [|3|,|4|.|lT 



This structure is more general than in perturbation theory and is the essential ingredient 
that allows SSB to coexist with "triviality". Finally, the physical normalization condition 

sua 

rf 2yl-loop | ^ 2 

—£jjr—kn=d*n- m h ( 9 ) 

determines 



and leads to 



(10) 



ml = M\ 2 R (11) 

In the context of the standard model, where Vr is known from the Fermi constant to be 246 
GeV, this predicts a Higgs mass ~2.2 TeV (up to radiative corrections that are small if the 
top mass is below 200 GeV 10). 

As discussed in |8|,p!0|,|3|] , one expects Eq.(ll), obtained in the single-component theory, 
to be also valid in the 0(N)-continuous symmetry case. This observation originates in Ref. 
|j~5| which obtained the same effective potential for the radial field as in the one-component 
theory. This is extremely intuitive. The Goldstone-boson fields do not contribute non- 
trivially to the effective potential: they contribute only their zero-point energy, which is a 
constant since these are free, massless fields, according to "triviality". Thus, in the 0(2) 
case, one may take the diagram (Vj=//,0s) for the one-component theory and "rotate" it 
around the V e ff symmetry axis. This generates a three- dimensional diagram (V^, 0i, 02) 
where V e ff depends on the bare radial field, 

pb = \l<t>\ + <f>l (12) 

in exactly the same way as V e ff depends on (fts in the one-component theory; namely 
{u\p B ) = 3X oP l) 

T/l-loop/ N A ° 4 , UJ 4 (p B ) f UJ 2 (p B ) l\ 

V (p*)-^Pb + ^2^—-- 2 }. (13) 

This represents the classical potential plus the zero-point energy of the free shifted radial 
field. "Triviality" implies that there are no observable interaction effects, so this result 
should be exact |||4j]. By using eqs.(6,10), y 1-10013 can be re-expressed in the form 

V^(PB) = ^(\n P 4-h (14) 
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Differentiating Eq.(14), we obtain the bare "radial source" 

j{pb)= h^m = ^a (15) 

dp B ^ v b 

which we shall compare with the lattice results for J = J(pb)- 

The lattice simulation of the 0(2)-invariant theory is obtained from the action 

4 ^ r S^i^i^ + ae^-^Cx)) 2 E^ig^ + ae„) - $ 2 (x)) 2 
° ^ [ 2a2 + 2a 2 

+ ^r G ($ 2 (x) + $ 2 (x)) + ^($ 2 (a;) + $ 2 (a;)) 2 ] (16) 
and one couples $i and $2 to two constant external sources J\ and J2 through 

a 4 £ [Ji$i(x) + J 2 $ 2 (x)] (17) 

X 

By using J\ = J cos and J2 = J sin (9 it is straightforward to show that the bare radial field 

(01 = ($1)^2, 02 = ($2)ji,j a ) 



PB = V0i + 0i (is) 



does only depend on J, that is 



Pb = P b{J) (19) 



We started our Monte Carlo simulation on a 10 4 lattice by investigating first the (r , A ) 
correlation which corresponds to the classically scale-invariant case. Analytically, this cor- 
responds to determine r a from the zero-mass renormalization condition [|16| 



d 2 v f 



e// | 

„2 \pb=0 



(20) 



dp 2 B 

so that the theory does not contain any intrinsic scale in its symmetric phase ($1) = ($2) = 
0. However, on the lattice, due to the relatively large errors introduced from the direct use of 
Eq.(20), it is more convenient to define the massless theory as in for the one-component 
theory. There, we started from the general expression 0] 



J(</> B ) = a<l>%hi(<j>l/vl)+Pv%<j> B (l - (21) 

which is still consistent with "triviality" (corresponding to an effective potential given by 
the sum of a classical background and the zero point energy of a free field) but allows for 
an explicit scale-breaking term f3. Setting a = one obtains a good description of the data 
in the "extreme double well" limit (r a much more negative than r c , where r c corresponds 
to the onset of SSB) where SSB is a semi-classical phenomenon and the zero-point energy 
represents a small perturbation. Then, we started to increase r , at fixed A , toward the 
unknown value r c and examined the quality of the fit with (a, /3, vb) as free parameters. The 
value of r Q at which the quality of the 2- parameter fit {at, (3 = 0,v B ) becomes exactly the 
same as that of the more general 3- parameter case was used to define the massless case. At 
X — 1 the massless regime of the one-component theory was found at a value r Q = r s where 
r s a 2 ~ —0.45. Finally, the accurate weak-coupling relation between the bare mass and the 
euclidean cutoff (T^j (A# = 6A C ) 

.. A * a 2 = -^-A 2 (22) 



32tt 2 16tt 2 

was used to generate the massless theory at different values of A D from the previously de- 
termined value of r s at A = 1. It should be noted that our numerical coefficient relating 
ultraviolet cutoff and lattice spacing for quadratic divergences (Ao ~ 4.87 ) agrees well 
with an independent analysis of lattice data presented by Brahm jlTj which gives (in the 
range A G < 10 ) Aa = 4.893 ± 0.003. Also, ref. |17| predicts the massless regime to corre- 
spond to r s a 2 = —(0.224 ± 0.001) for A Q = 0.5 in the infinite-volume limit. This implies 
r s a 2 ~ —(0.448 ± 0.002) for A Q = 1, in excellent agreement with our result. 

In the 0(2) case, one can use the results of pl|L?| by modifying Eq.(22) for the two- 
component case. This simply introduces a combinatorial factor of 4/3 so that Eq.(22) 
becomes 



A (23) 



4vr 2- 

Hence, we expect the massless case to correspond to r s a 2 ~ —0.6 for A G = 1. This was 
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confirmed to good accuracy by using the above- described fitting procedure to Eq.(21) (now 
with (p B replaced by p#). Thus, the identification of the massless regime on the lattice does 
seem to obey the simple scaling laws (22-23) and is under theoretical control. 

In the analysis of the one-component theory it was found that Eq.(3) was poorly 
reproduced numerically at small values of J {a 3 \J\ ~ 0.01 or smaller). As a consequence, the 
values of 4>b (and any higher-order Green's functions) extracted from the direct computation 
at J = were not reliable. Similarly, in the 0(2) case we find that at small J the exact 
^-independence of ps (see Eq.(19) ), is poorly reproduced and our data processing becomes 
unreliable. We therefore consider a "safe" region of J-values, Ja? > 0.05, in which the 
spurious ^-dependence is less than ±3%. Fitting the data to Eq.(15) we can infer the values 
of vb and and compare with Eqs.(6,10). 

Our numerical values for ap#(J) in the massless case, obtained with the Metropolis 
algorithm on a 10 4 lattice, are reported in Table I for A„ =1.0, 1.5 and 2.0. For each A D 
the corresponding r Q is computed by using Eq.(23) and our numerical input r s a 2 = —0.6 for 
A D =1 . The errors in Table I are essentially determined from the observed spurious variation 
of pB in the range < 9 < 2tt. As discussed above, this is a numerical artifact and should 
be considered a systematic effect of the Monte Carlo lattice simulation. It is reproduced 
with three random number generators consistent with the Kolmogorov-Smirnov test at the 
level O(10~ 4 ). At low J this systematic effect completely dominates the error; the statistical 
errors, after 30,000 iterations, are 4-5 times smaller [|18| . 



Table I also reports the vb and values obtained from the two-parameter fits to the 
data using Eq.(15). The resulting Z^ values agree well with the one-loop prediction in 
Eq.(10). To perform a more stringent test of the one-loop potential we next constrain to 
its one-loop value in Eq.(10) and make a precise determination of avs from a one-parameter 
fit to Eq.(15). This allows a meaningful comparison with the lattice version of Eq.(6) 

(^) 1 - loop = ^=exp(-g) (24) 

where we have identified the euclidean ultraviolet cutoff A — > yi being an a priori 
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unknown coefficient. By replacing Eq.(24) into Eq.(15) we determined yi from the one- 
parameter fit to the data at A = 1.5 and r Q a 2 = —0.9 to be yi = 2.44 ± 0.03 {-j^j = jf)- 
(As discussed in 0, one does not expect precisely the same numerical coefficient to govern 
the relation between euclidean cutoff and lattice spacing for both quadratic and logarithmic 
divergences, see below). In Table II we show the results of the one-parameter fits to the data 
at A Q =1.0 and 2.0 and the comparison with Eq.(24) for = 2.44 ± 0.03. It is apparent 
from Table II that the one-loop potential well reproduces the lattice data, as previously 
discovered in the one-component case ||. 

The value of yi obtained from the lattice simulation of the 0(2) massless A$ 4 theory is 
~17% larger than the value yi = 2.07±0.01 obtained in the one-component case @. This is 
due to our choice of fixing the value of the bare mass r s a 2 = —0.6 at A D = 1 with the simple 
combinatorial factor 4/3 discussed above and ignoring finite size corrections to Eq.(23), see 



| 17j . Also, in the 0(2) case, errors (as estimated from the spurious ^-dependence of ps) are 
larger than in the one-component theory and, moreover, the massless regime appears not 
so sharply identified on the basis of the (J, ps) correlation. Indeed, it is found in a narrow 
range of r a - values near r Q a 2 = —0.6 for A D = 1. Choosing instead r s a 2 ~ —0.58 would give 
y L ~ 2.07. 

The agreement between one-loop predictions and Monte Carlo data is not a trivial test 
of perturbation theory but, rather, represents a non perturbative test of "triviality". If 
perturbation theory were valid then the data should agree at least as well, if not better, 
with the leading-log formula based on the perturbative /3-function 



5A l irx L 



1 + TTf hi 



47T 2 apg 

(xll denoting an adjustable parameter). However, when we fit the A Q = 1.0, 1.5 and 2.0 data 
to this formula we find, respectively, for 17 degrees of freedom, (x 2 )ll =9, 44 and 133 (to 
compare with the values \ 2 =0.8, 3.3 and 7.3 obtained from the 1-loop one-parameter fits 
with Eq.(15) when Z<p is constrained to its value in Eq.(10)). Note that we are in a region 
where the "A log" term is not small; it is of order unity. Thus, the good agreement between 
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the data and the one-loop formula is not because the higher-order corrections, expected 
by perturbation theory, are negligibly small; it is because they are absent. Without the 
"triviality" argument, this would be an incomprehensible miracle. 

In conclusion, the good agreement between lattice simulation and eqs. (15,10,23,24) pro- 
vides definite evidence that the dependence of the effective potential on the radial field in the 
continuous symmetry case is completely consistent with our expectations |||4[]. This con- 
firms that Eq.(ll), up to small radiative corrections due to the gauge and Yukawa couplings, 
controls the relation between the Higgs mass and the Fermi constant in the standard model 
if SSB is generated through "dimensional transmutation" from a classically scale-invariant 
A$ 4 theory. The "triviality" structure we have checked, in which rrih is not proportional 
to "Afl" (which vanishes), implies that the Higgs, despite of its rather large mass, is only 
weakly interacting and would be free if the gauge and Yukawa couplings would be turned 
off. 
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TABLES 

TABLE I. The values of clpb(J) for the massless case are reported as discussed in the text. At 
the various values of A c and r we also show the results of the 2-parameter fits with Eq.(15) and 
the one loop prediction (10). 



Ja 3 


A D = 1.0 r Q a 2 = -0.6 


A = 1.5 r a 2 = -0.9 


A = 2.0 r Q a 2 = -1.2 


0.050 


0.4110 ±0.0130 


0.3850 ±0.0116 


0.3753 ±0.0091 


0.075 


0.4686 ± 0.0086 


0.4337 ±0.0086 


0.4176 ±0.0069 


0.100 


0.5133 ±0.0084 


0.4724 ±0.0065 


0.4517 ±0.0055 


0.125 


0.5495 ± 0.0063 


0.5048 ±0.0065 


0.4814 ±0.0054 


0.150 


0.5819 ±0.0060 


0.5332 ±0.0053 


0.5063 ±0.0042 


0.200 


0.6377 ±0.0047 


0.5812 ±0.0040 


0.5500 ±0.0037 


0.250 


0.6844 ± 0.0039 


0.6217 ±0.0034 


0.5875 ±0.0029 


0.300 


0.7256 ± 0.0030 


0.6571 ±0.0026 


0.6190 ±0.0029 


0.350 


0.7612 ±0.0029 


0.6892 ±0.0026 


0.6473 ±0.0026 


0.400 


0.7938 ± 0.0025 


0.7178 ±0.0025 


0.6731 ±0.0025 


0.450 


0.8246 ± 0.0024 


0.7435 ±0.0025 


0.6969 ±0.0023 


0.500 


0.8520 ± 0.0023 


0.7683 ±0.0021 


0.7193 ±0.0023 


0.550 


0.8785 ±0.0021 


0.7911 ±0.0021 


0.7399 ±0.0018 


0.600 


0.9029 ±0.0019 


0.8124 ±0.0020 


0.7593 ±0.0017 


0.650 


0.9261 ±0.0019 


0.8330 ±0.0019 


0.7778 ±0.0017 


0.700 


0.9481 ±0.0018 


0.8520 ±0.0019 


0.7953 ±0.0017 


0.750 


0.9693 ±0.0018 


0.8703 ±0.0017 


0.8119 ±0.0015 


0.800 


0.9893 ±0.0016 


0.8879 ±0.0015 


0.8274 ±0.0013 




= 25.0 ± 1.6 


= 16.4 ±0.7 


= 12.3 ± 0.4 




av B = (1.4lJ;j>) 10- 3 


av B = (1.8±0.6)10- 2 


av B = (5.6±0.8)10- 2 



l-loop = 2g 3 z l-loop = 1? 5 z l-loop 
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TABLE II. By using Eq.(15), we show the results of the 1-parameter fits for av B at A =1.0 
and 2.0 when is constrained to its one-loop value in Eq.(10). We also show the predictions from 
Eq.(24), (av B ) Th , for yi = 2.44 ± 0.03 as determined from the fit to the data at A G =1.5 . 



A D = 1.0 r Q a 2 = -0.6 


A 


= 2.0 r Q a 2 = -1.2 


= 26.32 = fixed 




i = 13.16 = fixed 


aVB = (7.09±0.11)10- 4 


av B 


= (3.79±0.03)10- 2 


X 2 _ 0.8 




X 2 _ 7.3 


d.o.f 17 




d.o.f 17 



(av B ) Th = (6.85 ± 0.09)10" 4 (av B ) Th = (3.89 ± 0.05)10" 2 



13 



